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J. D. Vergados 

Theoretical Physics Division, University of Ioannina, Ioannina, Gr 451 10, Greece 

(Dated: February 2, 2008) 

The recent WMAP data have confirmed that exotic dark matter together with the vacuum en¬ 
ergy (cosmological constant) dominate in the flat Universe. Thus the direct dark matter detection, 
consisting of detecting the recoiling nucleus, is central to particle physics and cosmology. Super- 
symmetry provides a natural dark matter candidate, the lightest supersymmetric particle (LSP). 
The relevant cross sections arise out of two mechanisms: i) The coherent mode, due to the scalar 
interaction and ii) The spin contribution arising from the axial current. In this paper we will focus 
on the spin contribution, which maybe important, especially for light targets. 

PACS numbers: 95.35.+d, 12.60.Jv 


INTRODUCTION 

The combined MAXIMA-1 QJ. BOOMERANG ;2j, DASI fjf], COBE/DMR Cosmic Microwave Background (CMB) 
observations Q, the recent WMAP data [S] and SDSS b] imply that the Universe is flat ?] and and that most of the 
matter in the Universe is dark, i.e. exotic. 

n b = 0.044 ± 0.04, n m = 0.27 ± 0.04, U A = 0.69 ± 0.08 

for baryonic matter , cold dark matter and dark energy respectively. An analysis of a combination of SDSS and 
WMAP data yields ( 6 j fl m ~ 0.30 ± 0.04(ler). Crudely speaking and easy to remember 

Ul b « 0.05, UIcdm ~ 0.30, U A « 0.65 

Since the non exotic component cannot exceed 40% of the CDM Q, there is room for exotic WIMP’s (Weakly 
Interacting Massive Particles). In fact the DAMA experiment (9] has claimed the observation of one signal in 
direct detection of a WIMP, which with better statistics has subsequently been interpreted as a modulation signal 
EH.These data, however, if they are due to the coherent process, are not consistent with other recent experiments, 
see e.g. EDELWEISS and CDMS 0 - It could still be interpreted as due to the spin cross section, but with a new 
interpretation of the extracted nucleon cross section. 

Supersymmetry naturally provides candidates for the dark matter constituents El] , .13“ 0- I n the most favored 
scenario of supersymmetry the LSP can be simply described as a Majorana fermion, a linear combination of the 
neutral components of the gauginos and higgsinos 0 ], 0 - 0 - In most calculations the neutralino is assumed to be 
primarily a gaugino, usually a bino. Models which predict a substantial fraction of higgsino lead to a relatively large 
spin induced cross section due to the Z-exchange. Such models have been less popular, since they tend to violate 
the relic abundance constraint. These fairly stringent constrains, however, apply only in the thermal production 
mechanism. Furthermore they do not affect the LSP density in our vicinity derived from the rotational curves. We 
thus feel free to explore the consequences of two recent models 0 , 0 , which are non-universal gaugino mass models 
and give rise to large hi ggsi no components. Sizable spin cross sections also arise in the context of other models, which 
have appeared recently |lS|. 0 -;20] (see also Ellis et al [2lJ for a recent update of such calculations). 

Knowledge of the spin induced nucleon cross section is very important since it may lead to transitions to excited 
states, which provide the attractive signature of detecting the de-excitation 7 rays in or without coincidence with the 
recoiling nucleus E3“l23- Furthermore it may dominate in very light systems like the 3 He, which offer attractive 
advantages 0 


THE ESSENTIAL THEORETICAL INGREDIENTS OF DIRECT DETECTION. 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational 
curves, it is essential to directly detect 0 , 0-0 such matter. Such a direct detection, among other things, may 
also unravel the nature of the constituents of dark matter. The possibility of such detection, however, depends on the 
nature of its constituents. Here we will assume that such a constituent is the lightest supersymmetric particle or LSP. 
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Since this particle is expected to be very massive, m x > 30 GeV , and extremely non relativistic with average kinetic 
energy T ~ 50ATeU(m x /100GeU), it can be directly detected |lj|-[25| mainly via the recoiling of a nucleus (A,Z) in 
elastic scattering. The event rate for such a process can be computed from the following ingredients: 

1. An effective Lagrangian at the elementary particle (quark) level obtained in the framework of supersymmetry 
as described , e.g., in Refs mil- 

2. A well defined procedure for transforming the amplitude obtained using the previous effective Lagrangian from 
the quark to the nucleon level, i.e. 

^e//=-^(xnVxOA^) (1) 

where 

Jx{Z) = N'y\[fy(Z) + ti(Z)r 3 + f° A (Z) l5 + f\(Z) l5 r 3 \N (2) 

The superscripts 0(1) refer to the isoscalar (isovector) components of the current. This form depends, of course, 
on the quark model for the nucleon. This step is not trivial, since the obtained results depend crucially on the 
content of the nucleon in quarks other than u and d. This is particularly true for the scalar couplings, which 
are proportional to the quark masses [IV] — J2f| as well as the isoscalar axial coupling. 

3. Nuclear matrix elements. l29l| —130| obtained with as reliable as possible many body nuclear wave functions. 
Fortunately in the most studied case of the scalar coupling the situation is quite simple, since then one needs 
only the nuclear form factor. Some progress has also been made in obtaining reliable static spin matrix elements 
and spin response functions f29j| — )30]. 

Since the obtained rates are very low, one wouldTike to be able to exploit the modulation of the event rates due to the 
earth’s revolution around the sun In order to accomplish this one adopts a folding procedure, i.e one 

has to assume some velocity distribution pil ST 0 liUzl-Hi for the LSP. In addition one would like to exploit 
the signatures expected to show up in directional experiments, by observing the nucleus in a certain direction. Since 
the sun is moving with relatively high velocity with respect to the center of the galaxy, one expects strong correlation 
of such observations with the motion of the sun |39l |40|. On top of this one expects to see a more interesting pattern 
of modulation as well. 

The calculation of this cross section has become pretty standard. One starts with representative input in the 
restricted SUSY parameter space as described in the literature for the scalar interaction |l£j,|4l| (see also Arnowitt 
and Dutta E^ b We will only outline here some features entering the spin contribution. The spin contribution comes 
mainly via the Z-exchange diagram, in which case the amplitude is proportional to Z 2 — Z 2 ( Z 3l Z 4 are the Higgsino 
components in the neutralino). Thus in order to get a substantial contribution the two higgsino components should 
be large and different from each other. Normally the allowed parameter space is constrained so that the neutralino 
(LSP) is primarily gaugino, to allow neutralino relic abundance in the allowed WMAP region mentioned above. Thus 
one cannot take advantage of the small Z mass to obtain large rates. Models with higgsino-like LSP are possible, but 
then, as we have mentioned, the LSP annihilation cross section gets enhanced and the relic abundance Q x h 2 gets 
below the allowed limit. It has recently been shown, however, that in the hyperbolic branch of the allowed parameter 
space 0 , 0 even with a higgsino like neutralino the WMAP relic abundance constraint can be respected. So, even 
though the issue may not be satisfactorily settled, we feel that it is worth exploiting the spin cross section in the 
direct neutralino detection, since, among other things, it may populate excited nuclear states, if they happen to be 
so low in energy that they become accessible to the low energy neutralinos [22j|-|23i|. 


RATES 


The differential non directional rate can be written as 


dRundir — 


p(0) TO 

m x ArriN 


da(u, w)|r<| 


( 3 ) 


where d<j(u 7 v) was given above, p( 0) = 0.3 GeV/cm 3 is the LSP density in our vicinity, m is the detector mass and 
m x is the LSP mass 
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The directional differential rate, in the direction e of the recoiling nucleus, is 0: 

p{ 0 ) m l y/u „ 

aKdir =- i - \v\v.e <a{v.e) — acr(u,v o( - -=. — v.e) 


m x ArriN 

where &(x) is the Heaviside function and: 


2tt 


fl r V 


V2 


du 


d<r(u, v) == 0/ ,, + T, S pi n Fii(u)} 


2 (p r bv) 2 

where u the energy transfer Q in dimensionless units given by 

Q 


u = —— , Qo = [m p Ab] 2 = 40 A 4//3 MeV 

Q o 


(4) 


(5) 


( 6 ) 


with b is the nuclear (harmonic oscillator) size parameter. F(u) is the nuclear form factor and Fn(u) is the spin 
response function associated with the isovector channel. 

The scalar cross section is given by: 


__ / dr \ 2 

s “ { ~iAP) ] p ’ x ° 


i + 4fr 

j s 


„S ( rr N 2 a 2 


(7) 


(since the heavy quarks dominate the isovector contribution is negligible). x 0 * s t ^ ie L^P-nucleon scalar cross 

section. The spin Cross section is given by: 

S 


1 


_ / r'r \ 2 spin > > 

— V / \ J ®r> v° ^spini S spin — 

Mp) P ’ X 3(1+ ^) 

J A 


f 0 

Ja \2 


S(u) 


( 8 ) 


S(u) » 5(0) = [(^-flo(O)) 2 + 2^n o (0)fii(0) + Hi(0)) 2 ] 
Ja Ja 


(9) 


The couplings f\ (/)() and the nuclear matrix elements fli(0) (fl o (0)) associated with the isovector (isoscalar) com¬ 
ponents are normalized so that, in the case of the proton at u = 0, they yield £ sp i n = 1. The proton cross section 
given by: 


a p P x° = 3cro l f°A + /Al 2 = 3cr 0 |a p | 2 , cr 0 = 2-(G F m p ) 2 = 0.77 x 10 38 cm 2 = 0.77 x 10 2 pb 


( 10 ) 


wile for the neutron (J sp ™ a = 3<ro| — f \\ 2 = 3cro|a n | 2 . Here a p and a n are the usual proton and neutron spin 
amplitudes (Sj. 

With these definitions in the proton neutron representation we get: 


Cpin = xS (0) 


( 11 ) 


5(0) = 


( —H„(0)) 2 + 2^Q„(0)Q P (0) + H 2 (0) 


( 12 ) 


where f2 p (0) and O„(0) are the proton and neutron components of the static spin nuclear matrix elements. In 
extracting limits on the nucleon cross sections from the data we will find it convenient to write: 


spin > 
v o Sspin — 


H 2 (0) 


' spin ^ 


^ n 


A 8 12 


(13) 


In Eq. m s the relative phase between the two amplitudes defined by 

yA d q Aq N 


q=u,d,s 


a,N = 


(14) 
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TABLE I: The static spin matrix elements for various nuclei. For 3 He see Moulin, Mayet and Santos l-46ll . For the other light 
nuclei the calculations are from DIVARI Hi For 73 Ge and 127 I the results presented are from Ressel et al [2?1 (*) and the 
Finish group et al Idsll (**). For 207 Pb they were obtained by the Ioannina team (+). Efflll . jT4 . 



3 He 

i9p 

29 Si 

23 Na 

73 Ge 

127 j* 

127j** 

207 p b + 

Do(0) 

1.244 

1.616 

0.455 

0.691 

1.075 

1.815 

1.220 

0.552 

Di(0) 

-1.527 

1.675 

-0.461 

0.588 

-1.003 

1.105 

1.230 

-0.480 

n P (o) 

-0.141 

1.646 

-0.003 0.640 

0.036 

1.460 

1.225 

0.036 

Dn(0) 

1.386 

-0.030 

0.459 

0.051 

1.040 

0.355 

-0.005 

0.516 

fJ'th 


2.91 

-0.50 

2.22 





l^exp 


2.62 

-0.56 

2.22 





li th (spin) 


0.91 

0.99 

0.57 





H-exp 










2s p Aq N = (N\q"/^ 5 q\N) (15) 

where is the nucleon spin and d q the relevant spin amplitudes at the quark level obtained in a given SUSY model. 
The isoscalar and the isovector axial current couplings at the nucleon level, f A , f\, are obtained from the corresponding 
ones given by the SUSY models at the quark level, f A {q), /A(?)> via renormalization coefficients g A , g\, i.e. f A = 
g () A f\{q), f\ = g\f\(q). The renormalization coefficients are given terms of A q defined above 0 , via the relations 

g° A = Au +Ad +As = 0.77 - 0.49 - 0.15 = 0.13 , g\ = Au-Ad= 1.26 

We see that, barring very unusual circumstances at the quark level, the isoscalar contribution is negligible. It is for 
this reason that one might prefer to work in the isospin basis. The static spin matrix elements are obtained in the 
context of a given nuclear model. Some such matrix elements of interest to the planned experiments are given in table 
HI The shown results are obtained from DIVARI [3(j, Ressel et al (*) [2jS|, the Finish group (**) [45| and the Ioannina 
team (+) [dof . [25$ . 

The spin ME are defined as follows: 

fip(0) = J —j— -< J J\ a z{p)\J J >- , U»( 0) = ^ 1 A J J\a z (n)\J J >- (16) 

where J is the total angular momentum of the nucleus and <j z = 2 S z . The spin operator is defined by S z (p) = 
i.e. a sum over all protons in the nucleus, and S z (n) = S z (i), i.e. a sum over all neutrons. 

Furthermore 

n 0 (o) = D p (0) + D n (0) , n 1 (o) = n p (o)-n„(o) (17) 


EXPRESSIONS FOR THE RATES 


To obtain the total rates one must fold with LSP velocity and integrate the above expressions over the energy 
transfer from Q m in determined by the detector energy cutoff to Q m ax determined by the maximum LSP velocity 
(escape velocity, put in by hand in the Maxwellian distribution), i.e. v esc = 2.84 vq with vq the velocity of the sun 
around the center of the galaxy(229 Km/s). 

Ignoring the motion of the Earth the total (non directional) rate is given by 

R = Rt(a, Qmin) ( 18 ) 


R = 


0) 


m 


(^h ) 2 VW)l*Lo A 2 + a s v p ™ Cs Pin ] 


m x o Am p g r {p) 


The SUSY parameters have been absorbed in R. The parameter t takes care of the nuclear form factor and the folding 
with LSP velocity distribution 000 (see table [HJ. It depends on Q m in, he. the energy transfer cutoff imposed 
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by the detector and a = [p r bv o\/ 2 ] 1 . 

In the present work we find it convenient to re-write it as: 


R = K 


Ccoh(A, )l r (AP) ) 


<? 1 


spin 

cr i 


(19) 


with cr i = 10 5 pb and 


K = 


P( 0) 


m 


100 GeV 


\J(v 2 ) cr i ~ 1.60 10 2 y 


P( 0) 


m 




0.3GeVcra -3 liG/ 280 kms~ 


and 


Ccoh(.A, fl r (^ 4 )) — 


100 GeV 


fi r (-*4) 

Hr(p) 


1 2 


^4 t C oh (- / 4) 5 Cspin (-^j f^r (-^4)) — 


lOOGeV 


/i r (^4) 


12 


MrO) 


L spm 


(A) 


A 


( 20 ) 


( 21 ) 


The parameters c co h(A, fi r (A)), c sp i n (A, /i r (A)), which give the relative merit for the coherent and the spin contri¬ 
butions in the case of a nuclear target compared to those of the proton, are tabulated in table |TT| for energy cutoff 
Qmin = 0, 10 keV. 

Via Eq. {0 we can extract the nucleon cross section from the data. 

We distinguish the following cases: 


• If the nuclear contribution comes predominantly from protons (Hi = Hq = H p ), S(u) « fI 2 and 


C spin 


3 


( 22 ) 


So knowing the nuclear matrix element one can extract from the data the proton cross section. One such 
example with negligible neutron contribution is 19 F (see table [I]) 


• If the nuclear contribution comes predominantly from neutrons ((Ho = —Hi = H ra ), one can similarly extract 
the neutron cross section from the data. 


• In many cases, however, the nuclear structure is such that one can have contributions from both protons and 
neutrons. The situation is then complicated and will be discussed below (see exclusion plots below). 

We have seen, however, that in going from the quark to the nucleon level we encounter the renormalization 
factors g° A = 0.1 , g\ = 1.2. Thus the isoscalar contribution is suppressed and 


5(0) * Of , (spin 


H? 

T 


Then the proton and the neutron spin cross sections are the same. 


Neglecting the isoscalar contribution and using H 2 = 1.22 and fl 2 = 2.8 for 127 I and 19 F respectively the extracted 
nucleon cross sections satisfy: 


spin 

a p,x° 

a L° 


C-coh (-^4? Mr (-^4) ) ^ 

Cspin (A,p r (A))\ H 2 


xlO 4 (A = 127) , « xlO 2 (A = 19) 


(23) 


It is for this reason that the limit on the spin proton cross section extracted from both targets is much poorer. 

The quantity c sp i n (A, g r (A)) ( sp i n , when the isoscalar contribution is neglected and employing H 2 = 1.22 (2.81) for 
12 '/ ( 19 F), is shown in Fig^I n the case of the spin induced cross section the light nucleus 19 F has certainly an 
advantage over the heavier nucleus 12 'I (see Fig. Q. For the coherent process, however, the light nucleus is no match 
(see Table Hill. 

If the effects of the motion of the Earth around the sun are included, the total non directional rate is given by 


R — A. [c C oh (A, (v4) )fJp ^.o(l T /i(tt, Q min')cosoP) j 


(24) 


and an analogous one for the spin contribution, h is the modulation amplitude and a is the phase of the Earth, which 
is zero around June 2nd. We are not going to elaborate further on this point since amplitude h depends only on the 
LSP mass and is independent of the other SUSY parameters. 
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TABLE II: The factors cl9 = c co h( 19,/x r (19)), sl9 = c sp j„(19, ^(19)), cl9 = c co h(73, /x r (73)), s73 = c sp i„( 73,/i r (73)) and 
c!27 = c co h(127, /rT.(127)), s!27 = c sp in(127, /i r (127)) for two values of Qmin- Also given are the factors s3 = c sp i n {3, /i r (3)) for 


Qmin 


m x (GeV) 

keV 


20 

30 

40 

50 

60 

80 

100 

200 

0 

t(3,s) 

2.3 

2.3 

2.3 

2.3 

2.3 

2.3 

2.3 

2.3 

0 

s3 

29.1 

20.6 

15.9 

12.9 

10.9 

8.3 

6.7 

3.4 

0 

t(19,c) 

1.153 

1.145 

1.138 

1.134 

1.130 

1.124 

1.121 

1.112 

0 

t(19,s) 

1.132 

1.117 

1.105 

1.096 

1.089 

1.079 

1.072 

1.056 

0 

cl9 

11465 

10478 

9423 

8499 

7702 

6451 

5539 

3212 

0 

sl9 

31.2 

28.3 

25.4 

22.8 

20.6 

17.2 

14.6 

8.4 

0 

t(73,c) 

2.238 

2.166 

2.094 

2.028 

1.967 

1.865 

1.785 

1.559 

0 

t(73,s) 

2.270 

2.223 

2.175 

2.129 

2.086 

2.012 

1.952 

1.771 

0 

c73 

225512 

261081 

278461 

284755 

284569 

274313 

258838 

186743 

0 

s73 

41.6 

47.7 

50.3 

50.9 

50.4 

47.7 

44.4 

30.8 

0 

t(127,c) 

0.984 

0.844 

0.721 

0.621 

0.542 

0.430 

0.358 

0.213 

0 

t(127,s) 

0.948 

0.796 

0.671 

0.574 

0.501 

0.401 

0.340 

0.220 

0 

cl27 

205674 

224676 

222547 

211216 

196895 

168585 

145173 

82424 

0 

sl27 

12.3 

13.1 

12.8 

12.1 

11.3 

9.7 

8.5 

5.3 

10 

t(19,c) 

0.352 

0.511 

0.592 

0.639 

0.667 

0.710 

0.720 

0.773 

10 

t(19,s) 

0.340 

0.489 

0.563 

0.606 

0.631 

0.669 

0.676 

0.720 

10 

cl9 

3500 

4676 

4902 

4789 

4546 

4075 

3557 

2233 

10 

sl9 

9.3 

12.4 

12.9 

12.6 

11.9 

10.6 

9.3 

5.8 

10 

t(73,c) 

0 

0.020 

0.119 

0.246 

0.363 

0.539 

0.651 

0.847 

10 

t(73,s) 

0 

0.0175 

0.105 

0.213 

0.311 

0.453 

0.539 

0.677 

10 

c73) 

0 

2313 

15295 

32947 

49559 

73463 

86339 

89290 

10 

s73 

0 

0.39 

2.5 

5.3 

7.9 

11.6 

13.4 

13.4 

10 

t(127,c) 

0.000 

0.156 

0.205 

0.222 

0.216 

0.191 

0.175 

0.109 

10 

t(127,s) 

0.000 

0.135 

0.177 

0.192 

0.190 

0.174 

0.165 

0.121 

10 

cl27 

0 

41528 

63276 

75507 

78468 

74883 

70964 

42180 

10 

sl27 

0 . 

2.2 

3.4 

4.0 

4.3 

4.2 

4.1 

2.9 



FIG. 1: On the left the quantity c sp in(A, fj, r {A)) (, ap in for the A=19 system is shown for two cut off values Qmin = 0, continuous 
curve, and Qmin = 10 keV, dotted curve. On the right the same quantity is shown for the A=127 system. The advantages of 
the lighter target, especially for light LSP, are obvious. 


SOME SUSY INPUT 


The most important input towards computing the event rate is the spin nucleon cross section. For orientation 
purposes will only incorporate the results of two recent calculations: 

• Ellis et al ]19j have given such cross sections in the allowed SUSY parameter space taking into account all 
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FIG. 2: On the left the proton cross section in pb is shown as a function of the LSP mass (in GeV) obtained by Ellis et al 
E3.il (logarithmic scale on both axes). On the right we show the results obtained by Hisano et al |2Q1. 



FIG. 3: The event rate per kg-y for the A=19 system on the left and for the A=127 on the right using the proton cross section 
of Ellis et al 23.0- Otherwise the notation is the same as in Fig. Q 


available experimental constraints and provided an update of such calculations l2ll . They have also calculated 
l48ll the spin proton cross section as a function of the LSP mass. Excluding a few data points we present their 
results in Fig. 0 


• In a recent paper Hisano, Matsumoto, Noijiri and Saito (2c|have obtained the proton cross section for Wino and 
Higgsino like LSP at the one loop level for relatively high LSP mass, where the tree contribution is negligible. 
Their results for the spin cross section are shown in the same figure, Fig. [21 

The event rate obtained with the proton cross section of Ellis et al is shown in Fig. |oJ while for the proton 

cross section of Hisano et al [2(| it is shown in Fig. 0J 


BOUNDS ON THE SCALAR PROTON CROSS SECTION 

Before proceeding with the analysis of the spin contribution we would like to discuss the limits on the scalar proton 
cross section. In what follows we will employ for all targets the limit of CDMS II for the Ge target |4^|, i.e. < 2.3 
events for an exposure of 52.5 Kg-d with a threshold of 10 keV. This event rate is similar to that for other systems 
(50| . The thus obtained limits are exhibited in Fiers I5l(il 


RESULTS-EXCLUSION PLOTS IN THE a p ,a n AND a p ,a n PLANES 

From the data one can extract a restricted region in the a p ,a n plane, which depends on the event rate and the 
LSP mass. Some such exclusion plots have already appeared j50|-|f>l| . One can plot the constraint imposed on the 




















FIG. 4: The event rate per kg-y for the A=19 system on the left and for the A=127 on the right using the proton cross section 
of Hisano et al (2C| . Otherwise the notation is the same as in Fig. 0 




GeV 


FIG. 5: The limits on the scalar proton cross section for A= 127 on the left and A= 73 on the right as functions of m x . The 
continuous (dashed) curves correspond to Qmin = 0 (10) keV respectively. Note that the advantage of the larger nuclear mass 
number of the A= 127 system is counterbalanced by the favorable form factor dependence of the A= 73 system. 


quantities \a p + fpOnl and|y4r^ + 


1 2 derived from the experimental limits via relations: 


fir 


| A ^ \J(Jji)C |“ A &bound{.A) v{ll T.^-,j4) . O'bound(A') — 


R 3 10 ~ 5 pb 


K SI 2 r 100 

^ ^ L p c spin 


( A. , Hr {A)) 


r(m x ,A) = 


cl° p ? n (A,Hr(A)) 


Cspin (A, Hr (A)) 


(25) 


where <5 is the phase difference between the two amplitudes and Cg°° n (A, /u r (A)) is the value of c sp i„(A, /z r (A)) evaluated 
for the LSP mass of 100 GeV. Furthermore 


| &p q <3.n| O'bound(A) [r(77r^- , a!)] ^ , &bound(^A) — 

i Lp 


&bound(A^ 


3cr 0 


1/2 


(26) 


The constraints will be obtained using the functions c^ n (A, Hr (A)), obtained without energy cut off , Q min = 0, 
even though the experiments have energy cut offs greater than zero. Furthermore even though we know of no model 




GeV 


FIG. 6: The same as in Fig. 0for A= 19 on the left and A= 3 on the right. As expected in the case of the coherent process 
the light systems are no match against the heavier ones (compare this figure with Fig. 0) 
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FIG. 7: The coefficients r(m x , 127) on the left and r(m x , 73) on the right as functions of m x . Note that due to the adopted 
normalization, for sufficiently large m x there is little difference between zero threshold (continuous curve) and a threshold of 
10 keV (dashed curve). For small LSP mass the results with a threshold of Qmin = 10 keV are not very reliable, since the 
relevant rate is very small. 




FIG. 8: The same as in figQfor the coefficients r(m x , 19) on the left and r(m x ,3). Only zero threshold was considered for 

3 He. 


such that e xS is complex, for completeness we will examine below this case as well. Such plots depend on the relative 
magnitude of the spin matrix elements. They will be given in units of the A-dependent quantity <Jbound(A) for the 
nucleon cross sections and the dimensionless quantity abound for the amplitudes respectively. The quantities r(m x , A), 
obtained from the data of table [HI are plotted for A = 127, 73, 19 and A = 3 in Figs. [7| and |8j Before we 
proceed further we should mention that, if both protons and neutrons contribute, the standard exclusion plot, must 
be replaced by a sequence of plots, one for each LSP mass or via three dimensional plots. We found it is adequate to 
provide one such plot for a standard LSP mass, e.g. 100 GeV, and zero energy threshold. The interested reader can 
deduce the scale for any other case with the help of Figs 0 and [HI In what follows we will employ for all targets the 
limit discussed in the previous section. 

Spin matrix elements of the same sign and Q, p 2> fin 

The situation is exhibited in Figs IDIlfl in the interesting case of the A=127 system using the nuclear matrix elements 
of Ressel et al given in Tabled One can understand the asymmetry in the plot due to the fact that fl p is much larger 
than O ra . In other words if a p happens to be very small a large a n will be required to accommodate the data. In the 
example considered here, however, the extreme values differ only by 20% from the values on the axes, which arise if 
one assumes that one mechanism at a time (proton or neutron) dominates. 

Spin matrix elements of opposite sign and |fi p | |fi n | 

This situation occurs in the case of the target 19 F'. As we have already mentioned the corresponding spin matrix 
elements shown in dare quite reliable. The obtained results are shown in Figs ll2ITH 
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FIG. 9: The boundary in the a p , a n plane extracted from the data for the target 12 ' I is shown assuming that the amplitudes 
are relatively real. The scale depends on the event rate and the LSP mass. Shown here is the scale for m x = 100 GeV. Note 
that the allowed region is confined when the amplitudes are of the same sign (left plot), but they are not confined when the 
amplitudes are of opposite sign. The allowed space now is i) The small triangle and ii) The space between the two parallel lines 
and on the right of the line that intercepts them. We also indicate by a dot the point a p = —a n favored by the spin structure 
of the nucleon. The nuclear ME employed were those of Ressel and Dean (see table [IJ 




5.0 x 10' 3 pb 


FIG. 10: The same as in Fig. [§]for the a p ,a n plane. On the left the allowed region is that below the curve (the amplitudes 
are relatively real and have the same sign) . In the plots on the right the amplitudes are relatively real and of opposite sign. 
The allowed region is i) between the higher segment of the hyperbola and the straight line and ii) Between the straight line 
and the lower segment of the curve. The nuclear ME employed were those of Ressel and Dean (see table QJ 



FIG. 11: The same as in Fig. EH assuming that the amplitudes are not relatively real, but are characterized by a phase 
difference 5. The allowed space is now confined. The results shown for the thin solid, thick solid and dashed curves correspond 
to <5 = 7t/3, 7t/6 and n/2 respectively . 
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FIG. 12: The same as in Fig. |H1 for the target 19 F. The matrix elements employed are those of tablcQ] Since now the proton 
and neutron spin ME have opposite signs, the conclusions about the relative sign of the elementary amplitudes a p ,a n are 
opposite of those arrived in Fig. |^1 The limit extracted on the neutron amplitude is poor, since the neutron spin ME is tiny. 



FIG. 13: The same as in Fig. ^]for 19 F with remarks about the signs of the elementary amplitudes as in Fig. m The matrix 
elements employed are those of table Q] 


Spin matrix elements of the same sign and 

This situation occurs in the case of the target 73 Ge (see table HJ. The obtained results are shown in Figs II 5lT7l 
Now the most stringent limit is imposed on the neutron cross section. 


Spin matrix elements of opposite sign and |ff n | |fi p | 

This situation occurs in the case of the target 3 He. The quantities t and h are essentially independent of the LSP 
mass, since the LSP is expected to be much heavier than the nuclear mass. We find t = 2.33 and h = 0.045 as the 
modulation amplitude. The latter is a bit larger than in heavier nuclei. This means that c sp i n (3, ft r (3)) =3 1 = 7.00. 
Likewise the parameters <Jbound{3) , r(m x , 3) are also constants, <Jbound,(3 ) =, r(m x ,3) = 1. Thus the expected rate 
is expected a bit smaller. This nucleus, however, has definite experimental advantages S3- Furthermore the nuclear 
matrix elements can be calculated very reliably. The differential tate is also independent of the LSP mass. Thus it 
can simply be exhibited in the form: 

= RL(u) [1 + H{u) cos a] (27) 

du 

Where L(it), integrated from 0 to u max = 1 yields unity, and H(u ) is the ratio of the modulated by the unmodulated 
differential rate. The situation is shown in Fig. 1181 The constraints on the elementary amplitudes and the nucleon 
spin cross sections are shown in Figs llOl£Tl Again the most stringent limit is imposed on the neutron cross section. 


CONCLUSIONS 


In the present paper we have studied the contribution of the axial current to the direct detection of the SUSY dark 
matter focusing on the popular targets 127 I, 73 Ge, 19 F and 3 He. The nuclear structure of these targets seems to favor 
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FIG. 14: The same as in Fig. 1111 for 19 F with remarks about the signs of the elementary amplitudes as in Fig. 1121 Note that 
now that the spin matrix elements are of opposite sign the extracted limits on the nucleon cross section are factor of two larger 
than assuming only one amplitude. 




200 


FIG. 15: The same as in Fig. [HI for the target 73 Ge. The matrix elements employed are those of table 0 The proton and 
neutron spin ME have the same sign, but now Q, n ftp. Thus the limit on a n is more stringent than that on a p . 


one component, either neutron, |U„| |fl p |, or proton, |fi p | |fi n | (see tabled . The real question is the size and 
the relative importance of the elementary amplitudes a p and a n . These result from the combination of two factors: 

• The relevant amplitudes at the quark level. 

In the case Z-exchange the the isoscalar amplitude is zero. In the case of the s-quark exchange the relative 
importance of the two amplitudes depends on the details of the SUSY parameter space. A priori there is no 
reason to prefer one amplitude over the other. 

• Going from the quark to the nucleon level. 

The structure of the nucleon tends to favor the isovector component. Thus, barring unusual circumstances 
at the quark level, the amplitudes a p and a n have opposite sign. Unfortunately none of the nuclear systems 
considered satisfies the most favorable condition fl n = — Nonetheless in this case one will be able to simply 
extract the amplitudes a p and a n from the experimental data, if and when they become available. 

In the most general case the extraction of the elementary amplitudes from the data is not going to be trivial. One 
can only make exclusion plots in the two dimensional ( a p , a n ) and (a p , cr p ) planes. By comparing such plots involving 
various targets one may be able to extract both amplitudes. In this work we have drawn such plots, using as input 
the nuclear spin response functions ( static spin ME and the relevant form factors obtained in the context of the shell 
model) in conjunction with the experimental limit, taken to be -< 16 events per kg-y for all targets considered. 
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FIG. 16: The same as in Fig. ^]for 73 Ge with similar remarks about the nucleon cross sections as in Fig. ll5lfor the amplitudes. 
The matrix elements employed are those of table 0 
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